Collective excitations of Dirac electrons in a graphene layer with spin-orbit interaction 
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Coulomb screening and excitation spectra of electrons in a graphene layer with spin-orbit inter- 
action (SOI) is studied in the random phase approximation. The SOI opens a gap between the 
valence and conduction bands of the semi-metal Dirac system and reshapes the bottom and top 
of the bands. As a result, we have observed a dramatic change in the long-wavelength dielectric 
function of the system. An undamped plasmon mode emerges from the inter-band electron-hole 
excitation continuum edge and vanishes or merges with a Landau damped mode on the edge of 
the intra-band electron-hole continuum. The characteristic collective excitation of the Dirac gas is 
recovered in a system with a high carrier density. 
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Since the recent discovery of quantum Hall effect in 
graphene [H [3- Q ■ there has been an upsurge of inter- 
est in understanding its electronic properties. Graphene 
has a promising potential for nanoscale device applica- 
tions [4J and is also very interesting physically because 
of its unusual Dirac- Weyl type band-structure near the 
Fermi level 0, 13, 0, S la E3, EH- As a single layer of 
graphite or an unrolled single-walled nanotube, the en- 
ergy band structure of graphene is well known just as 
for the graphite or for a nanotube H, 13 ■ The past 
few years have also evidenced increased activities in the 
role of spin-orbit interaction (SOI) in nanostructures be- 
cause it introduces many unusual features in these sys- 
tems. The SOI is also found to exhibit interesting effects 
such as the spin Hall effect in graphene B but overall 
its effect in graphene is as yet, unclear 12]. Further, 
the electron-electron interaction is important to under- 
stand the behavior of electrons in graphene and has been 
studied by many authors in gr aphi te-based structures but 
without the SOI included Gj,[HEl[H. In tnis letter > 
we explore the effect of SOI on the electron-electron in- 
teraction and the characteristic excitations in a graphene 
layer. 

Graphene has a honeycomb lattice of carbon atoms 
with two sublattices. Its energy band can be calcu- 
lated by the tight-binding model [3 and an intrinsic 
graphene is a semimetal with the Fermi energy located at 
the inequivalent K and K' points at opposite corners of 
its hexagonal Brillouin zone |5j . In the effective mass ap- 
proximation, the Hamiltonian of electrons near the Fermi 
energy is expressed by a 8 x 8 matrix. Since the SOI due 
to the atomic potential only mixes the states correspond- 
ing to the two sublattices, we can reduce the matrix to 
four independent 2x2 blocks. In the representation of 
the two sublattices, the Hamiltonian of a spin-up elec- 
tron near the K point of the reciprocal space reads 0, |j| 



H = vp ■ <j + A so a z = 



A so 



-ihv\7 
-A sn 



(1) 



with er = (cr x ,a y ,a z ) the Pauli matrices in the pseu- 
dospin space of the two sublattices and V ± = d/dx ± 



id I By. Here A so is the strength of the spin-orbit in- 
teraction (SOI) and v = 10 6 m/s is the 'light' veloc- 
ity of the Dirac electron gas. The eigenfunctions are, 
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for the state \k, 



the conduction band of energy Et 



while '5, (r) < i 

k \ -p ^ h cos(a fc /2) 



Ahrl sin(a fc /2) 



VA2 Q + H 2 V 2 k 2 
for the state 



|fe, — ) in the valence band E k = — ^jA 2 so + h 2 v 2 k 2 with 



L.2 _i_ JU2 
x T *y 



tan^fc = k y /k x , tanctfe = hvk/A so , and k = 

For a bare Coulomb scattering of two electrons at states 
|fe,A) and \p, Ai) into states \k + q, A') and \p — q, A' x ) 
respectively, the interaction matrix elements are 



A,A',Ai,Ai _ A, A 



Ai.A' 



u k,p 



= 9k (q)Mq)9p ' 



(2) 



Here vq = e 2 /(2tQeiq) is the two-dimensional Coulomb 
interaction (in Fourier space) with the high-frequency di- 
electric constant = 1 11], g^ x (q) is the interaction 
vertex, and the index A = ± denotes the two bands. It 
is straightforward to show that the RPA dressed inter- 
action matrix elements have the same form as the bare 
interaction matrix elements, i.e. Eq. J3J) 16]. As a re- 
sult, the dielectric matrix is expressed as a unit matrix 
multiplied by a dielectric function 



e(q,uj) = 1 - u (<7)n(q,w) 
with the electron-hole propagator 

f[E X k+q ] f[E x ] 



n(q,w) = 4 



E 

A,A',fc 



\gl' x \qt 



hu - 



E k+ q 



iS 



(3) 



(4) 



The factor 4 comes from the degenerate two spins 
and two valleys at K and K'; and the vertex fac- 
tor reads \g^' X (q)\ 2 = [1 + AA' cos ctk +q cos a k + 
XX' sin ctk+q sin ak{k + q cos 9) /\k + q\\/2 with 9 being the 
angle between fe and q. Since the intra-band backward 
scattering at q = 2k and the inter-band vertical transi- 
tion at q = are not allowed under Coulomb interaction 
in the system, we have \g£ ~ A (0)| 2 = \g% A (2fe)| 2 = 0. 
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FIG. 1: Real (e r ) and imaginary (ej) parts of the dielectric 
function vs lu at T = (left panels) and T = 2 K (right panels) 
for an intrinsic graphene (Ef = 0) and the SOI strength 
A so = (solid), 0.08 meV (dashed), and 0.1 meV (dotted). 
The wave vector is q = 0.05 x 10 5 cm -1 . 



The collective excitation spectrum is obtained by finding 
the zeros of the real part of the dielectric function e r . 
For convenience we denote each zero as a plasmon mode 
which may differ from the convention used in other places 
where some Landau damped modes are not counted since 
they do not have poles for e -1 . 

First we explore the properties of an intrinsic graphene 
where the net electron density is equal to zero and 
E F = 0. In Fig. n (a) and (b), we show the real (e r ) 
and imaginary (ej) parts of the dielectric function at 
T = for different SOI strengths (HI A so = (solid), 
0.08 (dashed), and 0.1 meV (dotted) at a wave vector 
q = 0.05 x 10 5 cm -1 . In the absence of the SOI, i.e. for 
A so — 0, the zero temperature electron-hole propagator 
is given by f[o(q, cu) = — q 2 j\^\J v 2 q 2 — lu 2 ], which has 
been previously obtained via the renormalization group 
theory [13, [j]f. We have e r = 1 (for u) > vq) above 
the inter-band electron-hole continuum (EHC) edge and 
£j = (for lu < vq) below it since only the intcrband tran- 
sition is allowed for electrons. With an increasing A so , 
the peaks of e r and e», which a re located at the edge of 
the inter-band EHC, lu h = 2y / A 2 so + h 2 v 2 q 2 /4, shift to 
higher energies. At the same time, for lu > u>h, £r in- 
creases continuously with A so . In the intrinsic graphene 
there is no plasmon mode at the zero temperature. 

At a finite temperature, the intra-band transition is 
allowed and contributes to the electron-hole propagator 
of Eq. Q while the inter-band contribution is reduced 
from that at zero temperature because of the electronic 
occupation of the conduction band. For A so = where 
u>l = uh, a £r dip at the intra-band EHC edge u>l = vq 
in the high energy side is formed as shown by the solid 
curve in Fig. ^ ( c ) an d two roots of the equation e r = 
or two plasmon modes may appear, one at lu — vq while 
the other has a higher energy. In Fig.Q](d), we observe a 
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FIG. 2: Plasmon spectrum (thick curves) of an electron gas 
in an intrinsic graphene (Ef — 0) at T = 2 K with A so = 
in (a) and 0.08 meV in (b). Intra- (dark shaded) and inter- 
(light shaded) band single-particle continuums are also shown. 
The vertical dotted lines indicate the q values for which the 
dielectric function is shown in Fig. 0(c). 



finite value of ej for lu < lul due to the intra-band tran- 
sition and a reduced for lu > lul corresponding to the 
weakening of inter-band scattering as the temperature in- 
creases. As a result, the plasmon mode at vq is strongly 
damped while the other is weakly damped. With increas- 
ing A so , the e r peak due to inter-band transitions shifts 
with loh to a higher energy while the e r dip of intra-band 
transitions stays with lul- Since the e r peak has a lower 
energy than the e r dip initially at A so = 0, the peak and 
the dip will merge at first and split again. As a result, the 
e r curve is deformed in such a way that two extra zeros of 
e r or two new plasmon modes emerge at a finite A so . If 
A so is further increased the e r peak finally disappears as 
the gap between the e r dip and peak becomes wider and 
the inter-band interaction diminishes. Corresponding to 
the separation of the e r peak from its dip, the curve 
develops a gap between lol and lor [Fig-D(d)]. 

The uj versus q spectrum of the plasmon modes at 
T = 2 K is plotted by thick solid curves for A so = 
in Fig. El (a) and A so = 0.08 meV in Fig. El (b). We 
observe one strongly Landau-damped plasmon mode at 
lu = vq and a weakly damped one in the inter-band EHC 
with an approximate dispersion of lu cx ^fq near q = 
in the A so = case. A finite A so separates luh from 
lul and opens a gap between the intra- and inter-band 
EHC's. Near q = the former weakly damped plas- 
mon mode becomes undamped because it is now located 
in the gap. Due to a strong inter-band single-particle 
transition at the edge of the inter-band EHC, the dis- 
persion curve of this mode is squeezed to a lower energy 
when it approaches the inter-band EHC edge and then 
splits into three plasmon modes, i.e., two new plasmon 
modes emerge near the inter-band EHC edge. One of 



the modes remain undamped in the EHC gap (or the 
Ci gap) while the other two are located inside the inter- 
band EHC. The latter two modes merge and disappear 
at q near 0.07 x 10 5 cm -1 in Fig. [3(b). The undamped 
mode survives in larger wavevectors until it merges with 
the strongly damped mode at the intra-band EHC edge. 

As explained in Fig. ^ (c) , the emergence of undamped 
and extra plasmon modes is a direct result of the SOI in 
the graphene system, which opens a EHC gap and sep- 
arates the e r peak attributed to inter-band transitions 
from the e r dip attributed to intra-band transitions. The 
plasmon spectrum is a result of the interplay between the 
intra- and inter-band contributions to the dielectric func- 
tion and may change significantly as parameters of the 
system such as the temperature and the electron density 
vary as discussed in the following. 

By applying a gate voltage to a graphene layer, one can 
control the electron density and the Fermi energy in the 
system 0,0 and it is interesting to see the corresponding 
change in the excitations. Due to the symmetry of the 
band structure, systems having the same density of elec- 
trons or holes are equivalent. Here we consider a system 
with a net electron density and positive Fermi energy Ep. 
For a Dirac gas in a graphene without a SOI, the extra 
electrons in the conduction band reduce the inter-band 
scattering rate but enhances the intra-band scattering by 
increasing the length of the Fermi ring. At T = 0, a EHC 
gap of width 2v(kp — q) is opened above the intra-band 
EHC edge in the range < q < kp. For Ep > A so the e r 
has a dip in the higher energy side of uj — vq and a dimin- 
ished peak in the lower energy side similar to the case of 
E F = but for T > [Fig.[T](c)]. As for the correspond- 
ing ej, the peak in the higher energy side at uj = v q + + 
is flattened because the inter-band transition to the bot- 
tom of the conduction band is forbidden but a new peak 
appears in the other side at u> — vq — + due to the now 
allowed intra-band transition. As a result of the e r dip 
at uj = vq + + introduced by the intra-band transition, 
a Landau damped plasmon mode appears near uj = vq 
and another one above, with an approximate dispersion 
uj oc ^fq at the small q. As compared to the spectrum of 
an intrinsic graphene at a finite temperature [Fig. [2 (a)], 
here the latter mode is undamped near q = and has 
a flatter dispersion slope near the crossing of the disper- 
sion curve and the inter-band EHC edge. Instead of an 
emergence of the new plasmon modes when the plasmon 
dispersion curve exits the EHC gap of SOI the plasmon 
dispersion smoothly exits the EHC gap corresponding to 
a finite Fermi energy. 

In the presence of the SOI, the physical scenario 
changes from the A so = case for long-wavelength exci- 
tations if the Fermi energy is located near the SOI gap. 
In Fig. |21 we plot the dielectric function versus uj for dif- 
ferent Fermi energies at T = (left panels) and T = 2 
K (right panels) in a graphene with the SOI strength 
of A so = 0.08 meV. The dashed curves of Fig. ^ for 
the Ep = case is also plotted here for comparison. 
At T = 0, due to modification of the energy bands at 
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FIG. 3: Real (e r ) and imaginary (e^) parts of the dielectric 
function vs u at T = (left panels) and T = 2 K (right panels) 
in a graphene of A so = 0.08 meV with Ef = (dashed), 0.25 
meV (solid), and 0.4 meV (dotted), q = 0.05 x 10 5 cm" 1 . 



the bottom (or top), the intra-band EHC edge deviates 
from uj = vq to a lower energy yj A 2 so + h 2 v 2 (kp + q) 2 — 
a/ A 2 d + h 2 v 2 k' F for small q. In Fig. |3] (a) and (b), we 
observe a shift of the e r dip and of the edge of the lower 
Cj peak to an energy lower than vq for the Ep = 0.25 
meV curve (solid) . With the Fermi energy inside the con- 
duction band, the inter-band EHC edge at T = is ex- 
pr essed by ujjj w hen q > 2kp but shifts to a higher energy 
y/A 2 so + K 2 v 2 k\ + ^A 2 so + h 2 v 2 (k F -q) 2 for q < 2k F . 
As a result, the solid e r peak and the upper gap edge 
(kp = 0.036 cm -1 ) have higher energies than the dashed 
ones (kp = 0) in Fig. [3| (a) and (b). When the Fermi 
energy becomes higher, for example, as shown by the 
dotted curve of Ep — 0.4 meV in Fig. 0(a) and (b), the 
system reverts to the Dirac gas and the dielectric func- 
tion becomes similar to that in the A so = case. From 
the solid curve for e r in Fig. (a) we observe that the 
SOI shifts the inter-band e r peak to a higher energy and 
can change the weakly damped plasmon mode inside the 
inter-band EHC to an undamped mode in the EHC gap 
created by the SOI. 

At a finite temperature, however, the restriction to 
single-particle transitions by the Fermi energy is relaxed 
and both the intra- and inter-band EHC edges are the 
same as those for Ep — 0. In Fig. (c), where e r at 
T = 2 K is shown, the peaks and the dips of all the 
curves have the same energy and for different Fermi 
energies in Fig. (d) have the same gap. These charac- 
teristics of the single-particle excitation continuum edges 
are also shown in Fig. 0] by the light and dark shades for 
inter- and intra-band EHC's respectively. 

The effect of SOI on the zero-temperature plasmon 
spectrum of a graphene with a finite Fermi energy out- 
lined in Fig. [21 (a) is illustrated by its dispersion in Fig.Q] 
(a). At a wavevector q close to or higher than kp, strong 
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FIG. 4: The same as in Fig. H but with E F = 0.25 meV 
at T = (a) and T = 2 K (b). The solid curves are for 
A so = 0.08 meV and the dashed in (a) for A so = 0. The 
corresponding dielectric functions for the values of q indicated 
by the vertical dotted lines are given in Fig. [3] by the solid 
curves. 



inter-band transitions between states of similar pseu- 
dospin |fe, +) and |fe',— ) with k anti-parallel to k' are 
possible. The contribution of these transitions at a finite 
A so increases the real part of the dielectric function at 
high energies and moves one zero of e r to a lower energy. 
As a result, the corresponding plasmon dispersion curve 
is excluded from the inter-band EHC and confined in the 
EHC gap opened by the finite Fermi energy and the SOI. 
At a finite temperature T = 2K, the electrons spread 
to higher energies and the inter-band transition rate is 
reduced at large q but increased at a small q. In this 
case, the dispersion curve of the mode above vq remains 
almost intact except a splitting to three modes when it 



enters the inter-band EHC near q ~ 0.004 x 10 5 cm -1 as 
shown in Fig. 0] (b). Being different from the behavior 
exhibited in Fig. [3(b), here the two extra plasmon modes 
along u>h merge with each other soon after the splitting 
and emerge again at a higher q. 

In summary, we have derived the electron-hole propa- 
gator and the dynamically screened Coulomb interaction 
matrix of an electron system in graphene with spin-orbit 
interaction in the random phase approximation, using 
the standard and widely used technique for multicompo- 
nent systems. In a system of intrinsic graphene without 
the SOI, our result reduces to the analytical result of a 
Dirac electron gas previously obtained by the renormal- 
ization group theory. The spin-orbit interaction changes 
the Dirac gas semimetal to a narrow gap semiconduc- 
tor system. The carriers far from the gap remain as the 
Dirac type while those near the gap shows different char- 
acteristics. At a finite temperature, the SOI splits the 
inter-band single-particle continuum from the intra band 
one and an undamped plasmon mode exists in the gap of 
the single-particle excitation spectrum. For a gap of 0.16 
meV, this plasmon mode exists in a range of wave vec- 
tor of the order of 0.1 x 10 5 cm _1 , which could perhaps be 
observed in experiments. With net electrons or holes in a 
graphene, a plasmon mode above the intra-band single- 
particle continuum exists and the Coulomb interaction 
is screened at a high energy even at the zero tempera- 
ture. The SOI may push this mode to the single-particle 
continuum gap or split it into three modes. In contrast 
to the negligible effect of SOI to the dynamic Coulomb 
screening and the plasmon spectrum of a Fermi gas in a 
InGaAs quantum well ^?|]; the SOI may change signifi- 
cantly these properties of a Dirac gas in a graphene. 
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